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Valley-based electronics, known as valleytronics, is one of the keys to break through to a new stage of elec- 
tronics. The valley degree of freedom is ubiquitous in the honeycomb lattice system. The honeycomb lattice 
structure of silicon called silicene is an fascinating playground of valleytronics. We investigate topological 
phases of silicene by introducing different exchange fields on the A and B sites. There emerges a rich va- 
riety of topologically protected states each of which has a characteristic spin-valley structure. The single 
Dirac-cone semimetal is such a state that one gap is closed while the other three gaps are open, evading the 
Nielsen-Ninomiya fermion-doubling problem. We have newly discovered a hybrid topological insulator named 
the quantum-spin-quantum-anomalous Hall insulator, where the quantum anomalous Hall effect occurs at one 
valley and the quantum spin Hall effect occurs at the other valley. Along its phase boundary, single-valley 
semimetals emerge, where only one of the two valleys is gapless with degenerated spins. These semimetals 
are also topologically protected because they appear in the interface of different topological insulators. Such a 
spin- valley dependent physics will be observed by optical absorption or edge modes. 



Introduction: The intrinsic degrees of freedom of an elec- 
tron are its charge and spin, which lead to electronics and spin- 
tronics. The valley degree of freedom on honeycomb lattices 
is expected to provide us with the notion of valleytronics. Val- 
leytronics was originally proposed in graphene[QJ-|3]], where 
the states near the Fermi energy are 7r orbitals residing near 
the K and K' points at opposite corners of the hexagonal Bril- 
louin zone. The low-energy dynamics in the K and K' valleys 
is described by the Dirac theory. The valley excitations are 
protected by the suppression of intervalley scattering. How- 
ever it is hard to realize valleytronics in graphene since the 
gap is closed and since it is difficult to discriminate between 
the K and K' points experimentally. In this context, transition 
metal dichalcogenides(4j-l7|] become a new playground of val- 
leytronics, where a considerably large gap is open. In the most 
recent experimental progress, the identity of valleys manifests 
as valley- selective circular dichroism, leading to valley polar- 
ization with circularly polarized light, offering a possibility to 
a realization of valleytronics yHZfl. 

Recently, another honeycomb system of silicon named sil- 
icene has been ex peri mentally synthesized flsl— TlOh and the- 
oretically explored fllll - l 1411 . It consists of buckled sublat- 
tices made of A sites and B sites. Controlling the gap at 
the K and K' points independently by applying external 
fields such that electric fieldd, photo-irradiation II 1 311 and 
exchange field|14], we are able to generate a rich variety of 
topologically protected states in silicene, each of which has a 
characteristic spin- valley structure. 

Well-known topologically protected states are quantum 
spin Hall {OSH)|[l5|] and quantum anomalous Hall (QAH) 



insulators (1161 - 12 111 . They are characterized by the helical and 
chiral gapless edge modes, respectively, when nanoribbons 
are made ll22M24l1 . The QAH effect is the quantum Hall effect 
without Landau levels, while the QSH effect is the quantum 
Hall effect of spins rather than charges. 

In this paper, to search for a full control of the valley de- 
gree of freedom, we make the use of the silicene 's buckled 



structure and introduce different exchange fields Ma and Mb 
operating on the A and B sites. This could be realized by 
attaching ferromagnets to silicene from the front and back 
sides. We obtain rich phase diagrams as illustrated in Fig.l 
and Fig. 3. First of all, we are able to generate the spin- 
polarized QAH (SQAH) insulator together with single Dirac- 
cone (SDC) semimetals along its phase boundaries [Fig. 2]. 
The SDC semimetal is a remarkable state that has one mass- 
less Dirac cone and three massive Dirac cones, evading the 
Nielsen-Ninomiya fermion-doubling problem[25J. Second, 
a new finding is the quantum-spin-quantum-anomolous Hall 
(QSQAH) insulator. It is a new type of topological insulator 
such that, e.g., the QAH effect is realized at the K point while 
the QSH effect is realized at the K' point [Fig.4]. Third, an- 
other new finding is a single- valley (SV) semimetal such that, 
e.g., the gap is open (closed) at the K (K 1 ) point with spin 
degeneracy. It is different from the SDC state which has only 
one closed gap without spin degeneracy. These valley depen- 
dent band structures will be experimentally observed by spin- 
valley selective circular dichroism. In what follows we use 
notations s z =t>U tz = A, B, r] = K, K' in indices while 
s z = ±1, t z — = ±1 in equations for the spin, the 

sublattice pseudospin and the valley, respectively. 

Hamiltonian: Silicene is described by the four-band 
second-nearest-neighbor tight binding model JTH fl4il26 l. 
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where c\ a creates an electron with spin polarization a at site 
i, and (((i, j))) run over all the nearest (next-nearest- 
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neighbor) hopping sites. We explain each term, (i) The 
first term represents the usual nearest-neighbor hopping with 
the transfer energy t = 1.6eV. (ii) The second term repre- 
sents the spin-orbit (SO) coupling with Aso = 3.9meV, where 
a = (a Xj a yj cr z ) is the Pauli matrix of spin, with z/^ = +1 
if the next-nearest-neighboring hopping is anticlockwise and 
Uij = —1 if it is clockwise with respect to the positive z axis, 
(iii) The third term represents the second Rashba SO cou- 
pling with Ar2 = 0.7meV associated with the next-nearest- 
neighbor hopping term, where t l z = ±1 for i representing the 
A (B) site, and = d^ / | with the vector d^ connect- 
ing two sites i and j in the same sublattice. (iv) The fifth 
term[12] is the staggered sublattice potential term. Due to the 
buckled structure the two sublattice planes are separated by a 
distance, which we denote by 21 with i = 0.23A. It gener- 
ates a staggered sublattice potential oc 2£E Z between silicon 
atoms at A sites and B sites in electric field E z . (v) The forth 
term represents the sublattice dependent exchange magnetiza- 
tion with exchange field M t i (Ma or Mb)- Sublattice depen- 
dence may be induced by attaching two ferromagnets from 
both sides of silicene. 

The low-energy effective Hamiltonian derived from the 
tight binding model (Q]) is given by the Dirac theory around 
the Krj (K or K') point as 

H v = hv F (r]k x r x + k y Ty) + f]a z r z Xso ~ £E z r z 

) + Ma z + AM(j z T z , (2) 

where r a are the Pauli matrices of the sublattice pseudospin. 
We have defined the mean and staggered exchange couplings, 
M = \{M A + M B ) and AM = \(M A - M B ), v F = &at 
is the Fermi velocity, and a = 3. 86 A is the lattice constant. 

The coefficient of r z is the mass of Dirac fermions in the 
Hamiltonian (|2]), 

Al=Vs z *so + s z AM-£E z , (3) 

which may be positive, negative or zero. The band gap is 
given by 2 1 | . It is intriguing that we can control it by 
changing the electric field E z and the staggered exchange field 
AM. We may make a full advantage of the spin- valley depen- 
dent Dirac mass to materialize valley tronics in silicene. 

There are two types of topological phase transitions in 
silicene. One occurs when the Dirac mass becomes zero, 
= 0, while the other occurs when the tip of any band 
touches the Fermi energy. We discuss these two cases succes- 
sively. 

Phase diagram in (E Z ,AM) plane: We investigate the 
topological phase transition in the (E z , AM) plane by setting 
Ma = —Mb- The topological phase transition is governed 
solely by the Dirac mass A^ given by (O. It occurs when the 
sign of the mass changes. The topological quantum numbers 
are the Chern number C and spin-Chern number C s modulo 2. 
They are given by C = + Cf + Cf + Cf, and C s = 
\{C^ + Cf - Cf - Cf), where is the summation of 
the Berry curvature in the momentum space over all occupied 
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FIG. 1: (Color online) Phase diagram in the £E Z -AM plane. Heavy 
lines represent phase boundaries where the band gap closes. There 
are three types of topological insulators as indicated by QSH with 
(0, 1) and SQAH with (±1, |). There are two types of trivial band 
insulators as indicated by AF and CDW. There emerge the SDC 
semimetal and the SV semimetal in the phase boundary. A circle 
shows a point where the energy spectrum is shown in Figj2] 




FIG. 2: (Color online) The band structures of zigzagn silicene 
nanoribbons at the points indicated in the phase diagram [Fig[T|. The 
vertical axis is the energy in unit of t, and the horizontal axis is the 
momentum. We can clearly see the Dirac cones representing the en- 
ergy spectrum of the bulk. Lines connecting the two Dirac cones are 
edge modes. 

states of electrons with spin s z in the Dirac valley K^. It is 
straightforward to calculate , 

Cl = |sgn(A^J (4) 

as a function of the Dirac mass A^ . 

We present the topological phase diagram in the (E z , AM) 
plane. The phase boundaries are given by solving A^ = 0, 
which yields four heavy lines corresponding to s z =t4- and 
r] = K,K' in the phase diagram [FigHl. The topologi- 
cal numbers (C, C s ) are calculated with the use of ©. We 
illustrate the spin- valley structure of topologically protected 
phases in Fig 12 

First, there appear four types of insulators: (1) The spin- 
polarized QAH (SQAH) insulators with (C,C S ) = (±1, \). 
(2) The QSH insulator with (0, 1). (3) The trivial charge- 
density- wave-type (CDW) band insulator with (0, 0). (4) The 
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trivial antiferromagnetic-order-type (AF) band insulator with 
(0,0). We note that there are two types of trivial band in- 
sulators. The band gaps are different between the K and K' 
points in the AF insulator, while they are identical in the CDW 
insulator. 

Second, SDC metals appear in the three phase boundaries 
of the SQAH insulator. The SDC metal is originally found 
in silicene by applying photo-irradiation and electric field 
simultaneously! 13]. It is interesting that the SDC state is also 
obtainable without photo-irradiation. On the other hand, the 
S V semimetal appears at the point where the four topological 
insulators meets. They are topologically protected semimet- 
als, since they appear in the interface of different topologi- 
cal insulators each of which is topologically protected against 
small perturbations. 

We may apply an inhomogeneous electric field lfl2ll 
E z (x,y) or generate a domain wall in the antiferromagnet 
AM {x, y), which makes the Dirac mass inhomogeneous. For 
simplicity we assume the homogeneity in the x direction. The 
zero modes appear along the line determined by A] z (y) = 0, 
when A] z (y) changes the sign. We may set k x =constant due 
to the translational invariance along the x axis. We seek the 
zero-energy solution. The particle-hole symmetry guarantees 
the existence of zero-energy solutions satisfying the relation 
ipB = i£ipA with £ = ±1. Here, ipA is a two-component am- 
plitude with the up spin and down spin. Setting ip a (%,y) = 
e ik * x (j> A (y) 9 we obtain H 7] ifj A (x,y) = E r] ^ A (x,y), to- 
gether with a linear dispersion relation = f]^hv F k x . We 
can explicitly solve this as 
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FIG. 3: (Color online) Phase diagram in the Ma- Mb plane. Dotted 
lines represent the points where the band gap closes. Heavy lines rep- 
resent phase boundaries. There are five types of topological insula- 
tors as indicated by QAH with (±2, 0), QSQAH with (±1,1 /2) and 
QSH with (0, 1). There is one trivial band insulator as indicated by 
AF. There emerge the mQAH metal and the TM in the phase bound- 
ary. A circle shows a point where the energy spectrum is calculated 
and shown in Fig|4] 




where C is the normalization constant. The sign £ is deter- 
mined so as to make the wave function finite in the limit \ y\ —> 
00. This is a reminiscence of the Jackiw-Rebbi mode[28] pre- 
sented for the chiral mode. The difference is the presence of 
the spin and valley indices in the wave function. 

Phase diagram in (M A ,M B ) plane: We proceed to de- 
rive the topological phase diagram in the (Ma, Mb) plane 
by setting E z = 0, and make its physical interpretation. We 
present our result in Fig|3] where the phase boundaries are not 
determined by the condition A^ = 0. This condition gives 
two dotted lines in the phase diagram [Fig m. We show soon 
how to derive the phase boundaries. We have also calculated 
the band structure of a silicene nanoribbon with zigzag edges, 
which we give in FigHfor typical points in the phase diagram. 

Diagonalizing the Hamiltonian © we obtain four energy 
levels, 



El = ±V(aA R fc) 2 ±F ± (fc) 2 



(6) 



with F±(k) = r]s z M ± \J(hv F k) 2 ± (^A S o + AM) . When 
two energy levels coincide, the band gap becomes zero. This 
occurs at the K and K' points, where k± = 0. The band 
closes when AM = —r]Xso for the K v point. 



FIG. 4: (Color online) The band structure of a silicene nanoribbon at 
marked points in the phase diagram [Figj3]. The vertical axis is the 
energy in unit of t, and the horizontal axis is the momentum. Lines 
connecting the two Dirac cones are edge modes. 



We determine the phase boundary. It is determined by the 
condition that the tip of any band touches the Fermi energy. 
Silicene is a metal on the phase boundary. The energy spec- 
trum at k = is given by 



E 



rjs z t z X S o + Ms z ± AMs z t z 



(7) 



It vanishes when Ma = — Aso or Mb = Aso for the K point, 
and Ma = Aso or Mb = —Aso for the K' point. 

The topological numbers are C = C K ± C K and C s = 
Cf + Cf , where C = C\ + C£ and Q = \^{C\ - C%) 
with 



-t z = ^sgn (77*3 A so - 



M, 



(8) 



We illustrate the spin- valley structure of typical states in FiglU 
First, there appear four types of insulators: (1) The QSQAH 
insulators with (C, C s ) = (±1, 1/2), where the QAH effect 
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is realized at one valley in coexistence with the QSH effect 
at the other valley. It is intriguing that, e.g., the topological 
numbers (±1, 0) are assigned to the K point and (0, 1/2) to 
the K' point. (2) The trivial AF insulator with (0, 0). (3) The 
QSH insulator with (0, 1). (4) The QAH phases with (2, 0) 
for M > and (-2, 0) for M < 0. 

Second, at the boundary, the topological metal (TM) and 
the single-valley topological metal (SVTM) emerge. The 
band touches parabolically in them. Their emergencies are 
also protected topologically since they are sandwiched by dif- 
ferent topological insulators. 

Optical absorption: An interesting experiment to probe 
and manipulate the valley degree of freedom is optical 
absorption J2H7I |27|]. We briefly discuss that spin- valley de- 
pendent band structures we have found will be experimentally 
observable by spin- valley- selective circular dichroism. Cir- 
cular dichroism is a phenomena in which the response of the 
left- and right-handed circularly polarized light is different. 
To assess the optical selectivity of spin and valley by circu- 
larly polarized light, we compute the spin-valley dependent 
degree of circular polarization, between the top valence bands 
and bottom of conduction bands. If we neglect the Rashba 
terms (Ar2 = 0) we are able to obtain an analytic formula for 
the transitions near the K v point, 

\Pl(k)\ 2 = mlv 2 Al±r 1 I , (9) 

where P± (k) are the interband matrix elements of the left- 
and right-polarized radiation fields for spin s z at k, respec- 
tively, defined for a vertical transition from band u c to band 
u Y , as P± (k) = mo (u c (k) \ \ ^r^ \u w (k)). Especially we 
find the spin- valley dependent optical selection rule at k = 0, 

\Pl (0)| 2 = m\vl (1 i^sgn (A^)) 2 , (10) 

which will be detected experimentally. Such a spin-valley se- 
lective circular dichroism would lead to the eventual realiza- 
tion of spin- valley tronics. 

In this paper, making the use of the buckled structure, we 
have shown that we are able to make a full control of the spin 
and valley degrees of freedom in silicene. The emergence 
of a relevant state with a particular spin-valley structure is 
topologically protected. In exploring the phase diagram, we 
have found a hybrid topological insulator named Quantum- 
Spin-Quantum-Anomalous Hall (QSQAH) insulator, where 
two different topological insulators coexist: The QSH effect 
is realized at one valley while the QAH effect is realized at 
the other valley. The topological numbers are simply given 
by one half of the sum of those of the QSH and QAH insula- 



tors. Such a hybrid of two distinctive topological insulators is 
utterly a new state. 
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